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ABSTRACT 



. , Standard Model (SSM) including the constraints from universal soft supersymmetry break- 
On! 

O . ing and radiative breaking of the electroweak symmetry. The latter is enforced by using the 
one-loop Higgs effective potential which automatically gives the one-loop corrected Higgs 
^ ! boson masses. We perform an exhaustive search of the allowed five- dimensional parameter 
■ space and find that the neutralino relic abundance O^/iq depends most strongly on the 
ratio ^0 = 'n^o/'n^i/2- For ^ 1 the relic abundance is almost always much too large, 
whereas for ^ 1 the opposite occurs. For ~ 1 there are wide ranges of the remaining 
parameters for which O,^ ~ 1. We also determine that rriq ^ 250 GeV and m^- ^ 100 GeV 
are necessary in order to possibly achieve Q.^ ~ 1. These lower bounds are much weaker 
than the corresponding ones derived previously when radiative breaking was not enforced. 
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1. Introduction 



The fundamental observation by tlie COBE satellite of minute anisotropies in the 
cosmic microwave background radiation [|l|], has once again made evident the synergism 
between particle physics and cosmology [|2[, since the most compelling explanation for 
these observations appears to be found in inflationary theories of the early universe 
g]. Theoretically, the observed large-scale structure of the universe requires the presence 
of a dominant cold dark matter (CDM) component in the energy density 0. However, 
CDM may not be enough to produce the observed small-scale structure [^, and recent 
re-evaluations of this problem suggest the presence of a sub-dominant hot dark matter 
(HDM) component as well, in the form of ~ few eV neutrinos. 

Even though the CDM component could have various origins from the particle 
physics point of view it is customary to ascribe it to an (9(100 GeV) stable particle, i.e., 
the lightest neutralino (x) (and also the lightest supersymmetric particle) of the mini- 
mal supersymmetric Standard Model (MSSM) 0. The standard calculational procedure 
consists of searching through the vast parameter space of the model for 'cosmo logically 
favored' regions, that is, regions with near critical relic abundance of neutralinos, i.e., 
where ^ /iq, with = Px/ Po the fraction of the present total energy density (as- 

sumed to be the critical density) in the form of neutralinos, and ho the Hubble parameter 
in units of lOOkms"^ Mpc~^ (current observations indicate that 0.5 < ho < 1). A more 
conservative [i.e., cosmological-model-independent) outcome of this procedure is the iden- 
tification of regions of parameter space which are cosmologically disfavored, i.e., regions 
where ^x^o ^ ^o- 

Since depends on the annihilation cross section xx ~^ ^ill, the whole set of couplings 
and masses in the model need to be specified before Qx can be accurately computed. In 
the MSSM twenty-one parameters are needed to specify the relevant quantities. To make 
things tractable, several ad-hoc assumptions about the model parameters are usually made 
[|^,|10|, resulting in approximate results of some interest but of limited scope. Things 
simplify dramatically when, following well-established theoretical prejudices, the MSSM 
is embedded in a generic supergravity model with universal soft supersymmetry breaking 
TT],p!^,P!B||T^,P3| , since then the parameter count drops to just eight. The final step in 



this theoretical round-up is to enforce the requirement of radiative electroweak symmetry 
breaking, which cuts down the parameter count to five variables, introduces a correlation 
between the usual relic density variables (M2, p), and determines the Higgs boson masses 
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in terms of the basic model parameters. (Hereafter we refer to this unified model as 
the Supersymmetric Standard Model (SSM).) The latter two novel effects have important 
consequences in the determination of the allowed region of parameter space and in the value 
of ^x^o throughout this region. This ultimate calculation has been performed before for a 
'no-scale'-inspired supersymmetry breaking scenario (where mi/2 ^ ^q,A ^ 0) imposing 
the radiative breaking constraint at tree-level fl^ and one-loop [|TT| , although in both cases 
with the inaccurate approximation of tree-level Higgs boson masses.0 The purpose of this 
note is to present the results of this calculation in the context of supergravity models with 
universal soft supersymmetry breaking and radiative electroweak breaking in the one-loop 
approximation [jl^, including the ensuing one- loop corrected Higgs boson masses [lT9| , pO| . 

Besides the particle physics model used to calculate the annihilation cross section, 
various approximations have been used regarding the composition of the x state, the 
dominant annihilation channels, and the solution of the Boltzmann equation needed to 
determine O^/iq- See Refs. ||T3| , pT| , pr7| for a discussion of these matters. The actual relic 
density calculations in this paper have been carried out following the methods previously 



described in Ref. \13 



2. The favored regions 

The Supersymmetric Standard Model (SSM) is an SU{3) x SU{2) x [/(I) model with 
the minimal three generations and two Higgs doublets of matter representations, and which 
is assumed to unify into a larger gauge group at a unification mass of Mjj ^ 10^^ GeV. 
The parameter space of this model can be described in terms of three universal soft- 
supersymmetry breaking parameters: mo, A, the top-quark mass mt, and the ratio of 

Higgs vacuum expectation values tan (3 = V2/V1; the sign of the superpotential Higgs mixing 
term fi is also undetermined [ffBl. Several consistency and phenomenological constraints 



restrict the range of the model parameters []T^,|T^. For ease of comparison with earlier 
work we calculate O^^o P^^t it in the traditional relic density (M2,/i) plane for a 
given tan/3, where M2 = (a2/«c/)wi/2 ~ 0.83mi/2 ~ O.dOrrig is the SU{2) gaugino mass. 
This is accomplished by numerically inverting the relation fu, = ^{jrit). We also take 



^ Concurrently with our calculation, there has appeared a new calculation [[T^ imposing ra- 
diative breaking at tree-level but including the one-loop corrections to the Higgs boson masses, 
although in a more restrictive supersymmetry breaking scenario where A is determined, and also 
assuming \i,{Mu) = Xt{Mu). 
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tan /3 = 2, 8, 15, 30, = mo/m^/2 = 0.0, 0.2, 0.4, 0.6, 0.8, 1.0 and 2, 4, 6, 8; U = = 
-?o, 0, +^o; and vary < M2 < 300 GeV and < < 800 GeV so that rrig ^ 1 TeV and 
fine-tuning of tlie parameters occurs at two-orders-of-magnitude or less [|l^ . Tlie resulting 
set of figures gives a good sampling of the whole parameter space. In Fig. 1 we present a 
representative set of plots for ^^^q in the (M2, fx) plane for tan/3 = 2, 8, ^0 = 0.0, 0.4, 1.0, 

= 0, and both signs of fx. Note that since we have restricted M2 < 300 GeV, then 
is also restricted: for large M2, < Mi = ^tan"^ 9wM2 ~ ^ 150 GeV. This upper 
bound on M2 also eliminates the possibility of nearly pure higgsino x's.i 

The first evident feature of these figures is that the allowed parameter space is con- 
siderably constrained. For ^0 ~ 1 the radiative breaking constraints exclude the lower 
diagonal portion of the figures due to the determination of /x. For example, at tree-level 

/"Le = ^?/2 [Xl/2 + ^0^0 - i(Mz/mi/2)2] , (2.1) 

and clearly for ~ 1? ~ "^1/2 since Xi/2,Xo are 0{1) functions [|18|. This correlation 
between fj, and M2 is a distinctive feature of supergravity models with radiative electroweak 
breaking and (for ^0 ~ 1) eliminates a large portion of the otherwise allowed parameter 
space (for example, compare Fig. Ic with Fig. 2 in Ref. ^1^)- 

Another consequence of 

the radiative breaking constraints is the determination of the Higgs boson masses in terms 
of the model parameters. This implies that as one moves around the (M2, fj.) plane in Fig. 
1, the Higgs masses are varying continuously. In fact, with growing M2 {i.e., 'mi/2), the 
lightest Higgs boson (h) quickly approaches its upper bound (which is variable at one-loop 
but does not exceed 150 GeV for rUg < 1 TeV) and the other three Higgs bosons become 
quite massive. This means that: (i) the important Higgs-mediated s-channel annihilation 
into fermion pairs becomes suppressed due to massive propagators as M2 grows; and (ii) 
the XX ~^ hh channel is open throughout a fixed region in the (M2,/u) plane {i.e., where 
> rrih). These two effects suppress the annihilation rate and hence enhance the relic 
density relative to the usual calculations (with [jl2|,|l^,|4|,|l5| or without [|10| the supergravity 
mass relations) where the radiative breaking constraints are not imposed and all the Higgs 
masses are kept fixed (at moderate values) throughout the {M2, fx) plane (for example, 
compare Fig. Ic with Fig. 2 in Ref. |1^ 



For contours of and x-composition in the (M2,/i) plane see e.g., Fig. 1 in Ref. ^ 
^ We should point out that in Refs. [13,14| fx followed an opposite sign convention than in the 



present paper and Refs. |l6| , |20|Jlq] . 
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For lower values of the radiative breaking constraints apply as well but the 
phenomenological constraints bite more into the allowed region. The left edge of the 
allowed area is determined primarily by the LEP lower bound on the chargino mass 
(m^+ > 45 GeV) and the constraint rrit ^ 90 GeV. For ~ the allowed area is further 
suppressed ||T^ by demanding an electrically neutral lightest supersymmetric particle. 

Within the allowed area the value of fi^^o grows steadily with since larger values of 
J ~ m'^^^{ci +^o) suppress the important t- and w-channel sfermion-mediated XX ~^ ff 
if = annihilation channels. On the other hand, does not affect the relic density 
in any significant way since A mainly determines the degree of left-right stop mixing. 
However, affects the calculated value of fx (mainly for small fi) fl^, and therefore it 
affects the left edge of the allowed area; ^a > {^A < 0) moves the left edge to the left 
(right) . 

The effect of M2 ~ 0.83mi/2 ~ CSOm^ is to scale up all Higgs and sparticle masses 
and therefore one would expect ^x^o t° grow steadily with mi/2- Even though this is 
generally the case, the steady growth can be locally (in mi/2) depleted due to the presence 
of poles and thresholds of the annihilation cross section. The most prominent of these 
are the Z- and /i-poles (xx ^ Z,h ^ ff) which occur for ~ and ~ 

respectively, and the xX ~^ hh threshold for ~ nih- In Fig. 2 (top row) we show 
versus rrig for ^0 = I5 = 0, nit = 150 GeV, tan/3 = 2, 8, and both signs of fx (solid 
/X > 0, dashed ^ < 0). These 'rays' correspond to a path of increasing M2 and (and 
therefore m^) through the plots in Fig. Ic (with M2/|/i| ~ 0.5(0.7) for tan/3 = 2(8)) 
such that nit = 150 GeV. For tan (3 = 2 and fx > the Z-pole occurs for nig ^ 270 GeV 
and the h-pole does not occur (since > \mh always). For fi < Q the h- and Z-poles 
are close: ^ \fnh ~ 37 GeV for nig ^ 340 GeV and ~ for rrig ~ 380 GeV. 

For tan/3 = 8 the h- and Z-poles are even closer: ~ ^^z for rrig ^ 310 (370) GeV 
and ^ ~ 50 GeV for rrig ^ 320 (400) GeV for > (/U < 0). These features are 
evident in the figures. The xX ^ hh threshold is also noticeable as a drop in the rate 
of increase of n^hl. The de gree of effectiveness of this new annihilation channel depends 
on the x-composition. Pure bino x's do not couple to the h field. This is the case for 
the tan/3 = 2, > ray where at rrig ^ 440 GeV, ^ rrih ~ 70 GeV but only an 
almost imperceptible drop occurs. The other rays encounter the threshold when x is of 
mixed composition and therefore the expected drop is clearly observable. We remark that 
perhaps the most important effect of the one-loop corrections to the Higgs boson masses is 
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to shift the position of the poles and thresholds discussed above, relative to the tree-level 
case. 

Since mi/2 and determine to a large extent the magnitude of f^x^O; possible 
to scan the parameter space and determine the lowest values of these parameters which 
give fi^^o ~ ^0- have done this for h = 1/2 and find that it is necessary to have 
rriq ^ 250 GeV and m^- ^ 100 GeV to possibly achieve this. Lower values of m- ^ do not 
suppress the annihilation cross section sufficiently to give a large enough relic density. 
The bounds are substantially weaker than previously obtained (c./. rriq > 600 GeV and 
nil ^ GeV Ijl^) without imposing the radiative breaking constraints since (as discussed 
above) now the important Higgs-dependent annihilation channels also fade away with 
increasing mi/2- 

The effect of tan/3 on fi^^o ^^^^ X"Composition for fixed M2 and /x, which 

affects the annihilation rate directly. Heuristically we find that O^^o generally decreases 
for large values of tan jS. This variable also affects the size of the allowed area in the (M2, /u) 
plane. The allowed area first increases with tan/?, it peaks in size, and then it decreases 
for large tan (3 until it disappears completely [|18| . 



From Fig. 1 and the preceeding discussion it is clear that for ^0 ^ 1 it is quite possible 
to obtain O^^o ~ ^O' pluses, squares, diamonds for h = 1, 1/v^, 1/2. As ^0 decreases 
the fraction of the allowed area occupied by 'preferred' points decreases and for ^0 ^ 0.2 
it is no longer possible to obtain ~ for 0.5 < ho < 1. However, ^0 ^ 1 may still 

provide enough neutralino relic density for it to be the major component of the galactic 
halo. 



3. The excluded regions 

We have determined that for ^0 ~ 1 and rrig < 1 TeV there are no points in parameter 
space where fi^^o > 1 ^^d therefore one cannot conclusively exclude any of these points 
on cosmological grounds. Note however that if /i = l/-\/2 (1/2) then > 1/2 (1/4) will 

become disfavored. On the other hand, for ^0 ^ 2 there are always regions of parameter 
space where > 1 for almost all allowed values of M2. The fraction of the allowed 

area which becomes cosmo logically excluded increases with ^o- 

To explore this dependence in a quantitative way it is better to change variables from 
(M2, fx) to (M2, mt) for ^0 ~ 1- This is because for large ^0 values, the range of /i is much 
enlarged, i.e., fx ^ ^o"T'1/2 (see Eqn. (|2Tl])), whereas mt is always restricted to be below 
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~ 190 GeV [lT6| , p!8|] . We thus interchange /U with the basic parameter rrit. In Fig. 3 we show 
n^/ig the (M2, mt) plane for ^0 = 2, = 0, tan/3 = 2, 8, and both signs of Note the 
large number of cosmo logically excluded points (at least the stars). 

The dependences of ^x^o radiative breaking constraints, and on ^o,^a,^i/2j 

and tan (3 are qualitatively as discussed for the ^0 < 1 case above. The relic density 
increases steadily with ^q, except for values of M2 where poles and/or thresholds occur. In 
Fig. 2 (bottom row) we show O^/ig along a ray with mt = 150 GeV in Fig. 3. Compared 
to the top row, the bottom row has very similar features, except for the overall scaling of 
the vertical axis. The reason for this behaviour {i.e., the little change in m^) is that n 
does not change too much for ^0 = 1 ^ 2 and the contours in the (M2, fj) plane (see 
e.g., Fig. 1 in Ref. [1^]) tend to change little with /j, (for large enough /i); the change 
(and the corresponding shifts in m^) are larger for fj, < 0. Note that in contrast with 
the tan/3 = 2, ^ > 0, ^0 = 1 case, for ^0 = 2 the /i-pole does occur in this case for 
rrig ~ 150 GeV and ~ ^ruh ~ 27 GeV. Regarding the XX ~^ hh thresholds in this case 
(see Fig. 2 bottom row), it is interesting to note that for tan/3 = 2, > 0, the threshold 
effect, even though small it is nevertheless noticeable (relative to the corresponding top 
row case in Fig. 2). This is because the relative size of the drop in the annihilation cross 
section due to the new channel is much larger here (^0 = 2, small total cross section) than 
in the previous case (^0 = 1, large total cross section). 

To quantify the statement made above about larger values of ^0 giving larger cosmo- 
logically excluded fractions of the allowed region, we define fh = n{Q^hQ > h'^)/ntoti where 
n(OT^/iQ > /ig) is the number of points inside the allowed region which have f^^^o > ^0 
and ritot is the total number of points for a fixed discrete gridding of the (M2, mt) space 
(with Amt = AM2 = 5 GeV). In Fig. 4 we plot fh=i/2 versus ^0 for = — ^o^O, +^0? 
tan (3 = 2, 8, 15, 30, and both signs of fx. Due to the large extent of the computations to 
produce Fig. 4, these were performed in the tree-level approximation to radiative breaking 
and the Higgs boson masses (a full one-loop treatment would have required C(days) of 
NEC SX-3 supercomputer CPU time). However, in this particular instance we do not 
expect the full one-loop treatment to modify significantly the results in Fig. 4 (as veri- 
fied explicitly in some cases) because its main effect is to shift the positions of poles and 
thresholds, which basically do not affect the value of fh- As anticipated, as ^0 grows the 
fraction of excluded points grows, and would probably approach ^ 100% sooner were it 
not for the poles. In fact, from Fig. 2 one can see that increasing ^0 reduces the size of 
the interval in mg where O^/iq < ^0 i~ 0-25 in this case). Because these regions can only 
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shrink but not disappear completely, it is not possible to obtain an absolute upper bound 
on ^0 from cosmological considerations alone. However, large specific regions of parameter 
space can be ruled out in this way. 



4. Conclusions 

We have presented an exhaustive exploration of the parameter space of the SSM which 
due to the constraints of universal soft supersymmetry breaking and radiative breaking of 
the electroweak symmetry needs only five variables to be described. The latter constraint 
imposes novel relations between the traditional relic density variables and allows the de- 
termination of the Higgs boson masses in terms of the basic parameters of the model. 
These two new effects affect significantly the variation of O^/iq throughout the allowed pa- 
rameter space. Furthermore, our calculations include the state-of-the-art one-loop Higgs 
effective potential to enforce the radiative breaking constraints and obtain automatically 
the one-loop corrected Higgs boson masses. As such we believe this to be the ultimate 
calculational framework upon which more specific supergravity models would need to be 
investigated (see e.g., [p2[|). 

We find that O^/iq depends most strongly on the parameter = 5 as pre- 

viously observed [|14|. For ^ 1 the relic density is almost always (barring accidental 
depletion due to Z- and/or /i-poles) in conflict with current cosmological observations. For 
^0^1 the opposite occurs, although ^x^o ^till account for the dark matter in the 
galactic halo. For ~ 1 there is a wide range of the other parameters for which ~ 1- 
Besides these qualitative observations, one can draw more precise conclusions for specific 
points or regions of parameter space. It is also possible to obtain an approximate lower 
bound on the squark and slepton masses below which it is not possible to get Q^. ~ 1; 
we find rriq ^ 250 GeV and rrij ^ 100 GeV. These lower bounds are much weaker than 
those obtained previously without enforcing the radiative breaking constraints. Since 
on naturalness grounds |T^ we have limited rrig to be below 1 TeV, the value of 171^ is also 
restricted to ^150 GeV. These new bounds offer added hope that the supersymmetric 
spectrum may be discovered at the next generation of ee and pp machines. 
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Figure Captions 



Fig. 1. The neutralino relic density ^^Hq in the (M2, /u) plane for tan/3 = 2,8, = 0, 
both signs of /U, and (a) ^0 = 0, (b) ^0 = 0.4, (c) ^0 = 1-0. The various symbols 
denote the following ranges of n^hl: < 0.025 (dots), 0.025 - 0.10 (crosses), 
0.10 - 0.25 (diamonds), 0.25 - 0.50 (squares), 0.50 - 1.0 (pluses), > 1.0 (stars). 
For h = 1, 1/v^, 0.5 stars, stars+pluses,stars+pluses+squares are excluded on 
cosmo logical grounds. 

Fig. 2. The neutralino relic density ^x^o ^ ^ function of the gluino mass nig along rays 
with mt = 150 GeV in Fig. Ic (^0 = = 0) and Fig. 3 (^0 = 2,^^ = 0). 

(In Fig. Ic these rays correspond to M2/\fi\ ~ 0.5 (0.7) for tan/3 = 2 (8).) Solid 
(dashed) lines correspond to |U > (/U < 0). Note the various poles and thresholds 
of the neutralino annihilation cross section. 

Fig. 3. The neutralino relic density ^x^o {M2,mt) plane for ^0 = 2,^^ = 0, 

for tan/3 = 2,8 and both signs of fi. The symbols are as in Fig. 1. For 
h = 1,1/^^2,0.5 stars, stars+pluses,stars+pluses+squares are excluded on cos- 
mological grounds. 

Fig. 4. The fraction of the allowed parameter space in the (M2,mt) plane which is cos- 
mologically excluded for h = 1/2 {fh=i/2) versus ^0 for $,a = — ^0,0, +^0) both 
signs of iJ, and tan/3 = 2 (solid), 8 (dashed), 15 (dotdashed), 30 (dotted). 
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